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Abstract The Fermat quotient q p (a) := (a p 1 — 1 )/p, for prime p \ a, and the Wil- 
son quotient w p := {{p— 1)! + 1)//? are integers. If p\w p , then p is a Wilson prime. 
For odd p, Lerch proved that 1p{ a ) ~ w p)/p is a l so an integer; we call it the 

Lerch quotient i p . If p \ £ p we say p is a Lerch prime. A simple Bernoulli-number 
test for Lerch primes is proven. There are four Lerch primes 3, 103,839,2237 up 
to 3 x 10 6 ; we relate them to the known Wilson primes 5, 13,563. Generalizations 
are suggested. Next, if p is a non- Wilson prime, then q p (w p ) is an integer that we 
call the Fermat-Wilson quotient of p. The GCD of all q p {w p ) is shown to be 24. 
If p | q p (a), then p is a Wieferich prime base a; we give a survey of them. Taking 
a = w p , if p | q P {w p ) we say p is a Wieferich-non-Wilson prime. There are three up 
to 10 7 , namely, 2,3, 14771. Several open problems are discussed. 



1 Introduction 

By Fermat's little theorem and Wilson's theorem, if p is a prime and a is an integer 
not divisible by p, then the Fermat quotient of p base a, 

a p ~ x - 1 

q p {a):= , (1) 



and the Wilson quotient of p, 



w p := , (2) 



are integers. (See [25, pp. 16 and 19] and [26, pp. 216-217].) 

For example, the Fermat quotients of the prime p = 5 base a = 1,2,3,4 are 
q 5 {a) = 0,3, 16,51; the Fermat quotients of p = 3,5,7, 11, 13, 17, 19,23,29,31, . . . 
base a = 2 are 



1 
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1p{ 2 ) 



2 p-l 



1 



1,3,9,93,315,3855, 13797, 182361,9256395,34636833, . . . 



P 



[28, sequence A007663]; and the Wilson quotients of p = 2,3,5,7, 11, 13, 17, . . . are 

w p = 1 , 1 , 5, 103, 32989 1 , 36846277, 1230752346353, . . . 

[28, sequence A007619]. 

A prime p is called a Wilson prime [15, section A2], [25, p. 277] if p divides w p , 
that is, if the supercongruence 



holds. (A supercongruence is a congruence whose modulus is a prime power.) 

For p = 2, 3, 5, 7, 11, 13, we find that w p = 1, 1,0,5, 1,0 (mod p) (see [28, se- 
quence A002068]), and so the first two Wilson primes are 5 and 13. The third and 
largest known one is 563, uncovered by Goldberg [13] in 1953. Crandall, Dilcher, 
and Pomerance [4] reported in 1997 that there are no new Wilson primes up to 



Vandiverin 1955 famously said (as quoted by MacHale [21, p. 140]): 

It is not known if there are infinitely many Wilson primes. This question seems to be of such 
a character that if I should come to life any time after my death and some mathematician 
were to tell me that it had definitely been settled, I think I would immediately drop dead 
again. 

As analogs of Fermat quotients, Wilson quotients, and Wilson primes, we intro- 
duce Lerch quotients and Lerch primes in Section 2, and Fermat- Wilson quotients 
and Wieferich-non- Wilson primes in Section 3. We define them by combining Fer- 
mat and Wilson quotients in apparently new ways. 



2 Lerch quotients and Lerch primes 

In 1905 Lerch [20] proved a congruence relating the Fermat and Wilson quotients 
of an odd prime. 

Lerch's Formula. If a prime p is odd, then 



(p-1)! + 1=0 (mod p 2 ) 



5 x 10 s . 



p-l 

E #p( a ) = w p ( m ° d p) 



that is, 



p-l 



p-(p-l)! = (mod p 2 ). 



(3) 
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Proof. Replace a with ab in equation (1). Substituting a 1 '^ 1 = pq p (a) + 1 and 
b p = pq p (b) + 1, we deduce Eisenstein's logarithmic relation [10] 

q p {ab) = q p (a) + q p {b) (mod p) 

and Lerch's formula follows. For details, see [20] or [29]. □ 

Ribenboim [26, p. 218] explains the point of Lerch's formula this way: 

Since the Fermat quotient is somehow hard to compute, it is more natural to relate their 
sums, over all the residue classes, to quantities defined by p. 

Wilson quotients and Lerch's formula have been used (see [29]) to characterize 
solutions of the congruence 

1" +2" + ••• + *" = (k+l) n (mod/fc 2 ). 



2.1 Lerch quotients 

Lerch's formula allows us to introduce the Lerch quotient of an odd prime, by anal- 
ogy with the classical Fermat and Wilson quotients of any prime. 

Definition 1. The Lerch quotient of an odd prime p is the integer 

P P 2 

For instance, 

+ 3 + 16 + 51 -5 1 + 16 + 81 +256-5-24 

t*. = = = 13. 

5 5 25 

The Lerch quotients of p = 3,5,7, 11, 13, 17, 19,23,29, . . . are 

t p = 0, 13, 1356, 123229034,79417031713,97237045496594199, 
166710337513971577670,993090310179794898808058068, 
60995221345838813484944512721637147449, 

and for prime p < 62563 the only Lerch quotient £ p that is itself a prime number is 
£5 = 13 (see [28, Sequence A197630]). By contrast, the Wilson quotients w p of the 
primes p = 5,7, 11,29,773, 1321,2621 are themselves prime [15, Section A2], [28, 
Sequence A050299]. 
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We define Lerch primes by analogy with Wilson primes. 

Definition 2. An odd prime p is a Lerch prime if p divides i p , that is, if 

p-1 

Y,a p ~ l -p-{p-\)\ = Q (mod p 3 ). (4) 

a=\ 

Forp = 3, 5, 7, 11, 13, 17, 19, 23, 29,31, 37,41, 43, 47,53, 59, 61, 67,71, 73, 79, 83, 
89,97, 101, 103, . . . , we find that 

l p = 0,3,5,5,6,12,13,3,7,19,2,21,34,33,52,31,51,38,32,25,25,25, 
53,22,98,0,... (mod/?) 

[28, Sequence A197631], and so the first two Lerch primes are 3 and 103. 

We give a test for Lerch primes involving Bernoulli numbers. Ubiquitous in num- 
ber theory, analysis, and topology (see Dilcher [7]), they are rational numbers B„ 
defined implicitly for n > 1 by the symbolic recurrence relation 

(B + - B" +l =0. 

(Ribenboim [26, p. 218] says, "Treat B as an indeterminate and, after computing 
the polynomial in the left-hand side, replace B k by B^.") Thus for n — 1, we have 
(B+ l) 2 -B 2 = 2Bi + 1 = 0, and so B x = -1/2. Now with n = 2, we see that 
(B + 1 ) 3 - B 3 = 3B 2 + 3Bi + 1 = leads to B 2 = 1/6. In this way, we get 

11 15 

Bt, =Q.B 4 = ,B 5 =0,B 6 = —,B 7 =0.fi 8 = ,B 9 =0,Bio = —,.... 

30 42 30 66 

In 1937 (before the era of high-speed computers!) Emma Lehmer [19] showed 
that 5 and 13 are the only Wilson primes p < 21 1. To do this, she used her husband 
D. H. Lehmer's table of Bernoulli numbers up to B220, together with Glaisher's 
congruence [11] (see also [20]), which holds for any prime p: 

w p = B p -i + - - 1 (mod p) . (5) 

Here recall the definition 

- = (mod m) m I a, GCD(a,^) = 1. 

b 

For example, that 5 is a Wilson prime, but 7 is not, follows from the congruences 
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1 5 

w 5 = B 4 + -- 1 = -- =0 (mod 5), 
5 6 

w 7 =B 6 +i-l = -- ^0 (mod 7). 

Multiplying Glaisher's congruence by p and substituting pw^ = (p — 1) ! + 1 
yields £. Lehmer's test: A prime p is a Wilson prime if and only if 

pB p -\ =p-\ (mod p 2 ). 

We provide an analogous test for Lerch primes. 

Theorem 1 (Test For Lerch Primes). A prime p > 3 is a Lerch prime if and only if 

pB p -x = p+{p-l)\ (mod/? 3 ). (6) 

Proof. We first establish the following Criterion: an odd prime p is a Lerch prime 
if and only if 

(B + p) p = p 2 + p! (mod//). (7) 

To see this, recall the classical application of Bernoulli numbers called Faulhaber's 
formula (also known as Bernoulli's formula — Knuth [18] has insights on this): 

l" + 2" + - + (fc-l)"= ( + } - ■ (8) 

n + 1 

(See Conway and Guy [2, pp. 106-109] for a lucid proof.) Now set k = p and 
n = p - 1 in (8). It turns out that B p = (indeed, B 3 = B 5 = B n = Be, = ■ ■ ■ = 0; 
see [2, p. 109], [16, section 7.9]), and it follows that the congruences (4) and (7) are 
equivalent. This proves the Criterion. 

To prove the Test, note that for any odd positive integer p, the vanishing of B2it+i 
for k > 1 implies 

(B+py = pp+p-pi>- l B l+ £ (i k )p p ~ B »- 



k=\ 



The von Staudt-Clausen theorem [2, p. 109], [16, section 7.9], [25, p. 340] says in 
part that the denominator of B^k is the product of all primes q for which {q — 1 ) | 2k. 
(For instance, as (2—1) | 2 and (3 — 1) | 2, the denominator of Z?2 is 2 • 3, agreeing 
with Z?2 = 1 /6 ) Thus, if p is an odd prime, then on the right-hand side of (9) only 
fip_i has denominator divisible by p. From this we see, for p > 5, that p 4 divides 
the numerator of each term except p 2 B p _\. (For the k = (p — 3)/2 term, this uses 
p | ( p p 3 ) •) Therefore, the congruence 

(B + p)" = p 2 B p _ l (mod/) (10) 
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holds for all primes p > 3. Substituting (10) into Criterion (7) and dividing by p, 
we arrive at Test (6). □ 

As a bonus, (10) affords a proof of Glaisher's congruence. 

Corollary 1. The congruence (5) holds. Equivalently, if p is any prime, then 

P B p ^=p+{p-\)\ (modp 2 ). (11) 

Proof. To see the equivalence, substitute (2) into (5) and multiply by p. To prove (11), 
first verify it for p = 2 and 3. If p > 3, use (3), (8), and the fact that B p = to get 
{B + p) p = p 2 + pi (mod p 3 ). Then (10) and division by p yield (11). □ 

Notice that the congruences (6) and (11) are the same, except that in (6) the 
modulus is p 3 , while in (1 1) it is p 2 . However, one cannot prove Corollary 1 trivially 
(by reducing (6) modulo p 2 instead of /? 3 ), because (6) holds only for Lerch primes, 
whereas (11) holds for all primes. 



2.3 Computing Lerch primes: a surprising crossover 

Let us compare two methods of computing Lerch primes: Definition (4) and Test (6). 
Both require, essentially, computation modulo p 3 . The Test seems simpler, but on 
the other hand it requires computing B p -\ modulo p 2 . 
To find out which is faster, we used the code 

If [Mod[Sum[PowerMod[a,p-l,p~3] , {a,l,p-l}] - p - (p-1) ! , p"3] 
== 0, Print [p] ] 

in a Mathematica (version 7.0.0) program for (4), and we used the code 

If [Mod [Numerator [p*Mod [BernoulliB [p-1 ], p"2 ] - p - (p-1)!], p~3] 
== 0, Print [p] ] 

in a program for (6). Here Mod [a, m] gives a mod m, PowerMod [a, b, m] gives 
a b mod m (and is faster than Mod [ a"b, m] ), and BernoulliB [ k] gives B^. 

Table 1 shows the CPU time (on a MacBook Air computer with OS X 10.6 and 
2.13GHz Intel processor) for each program to decide whether p is a Lerch prime. 

Note the surprising crossover in the interval 10007 <p< 2001 1: before it, 
Test (6) is much faster than Definition (4), but after the interval the reverse is true. 
Notice also that for p > 10 4 the CPU times of (4) grow at about the same rate as p, 
while those of (6) balloon at more than double that rate. 

The programs for (4) and (6) searched up to 10 4 in about 47.3 and 0.6 seconds, 
respectively, and found the Lerch primes 3, 103,839, and 2237 (see [28, Sequence 
A197632]). There are no others up to 10 6 , by the program for (4), which consumed 
about 160 hours. (To run the program for (6) that far up was not feasible.) 

Marek Wolf, using a modification of (4), has computed that there are no Lerch 
primes in the intervals 1000003 < p < 44961 13 and 18816869 < p < 18977773, as 
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well as 32452867 < p < 32602373. His computation took six months of CPU time 
on a 64-bit AMD Opteron 2700 MHz processor at the cluster [17]. 



CPU time in seconds 



p 


Definition 


vs. 


Test 


5 


0.000052 


> 


0.000040 


11 


0.000069 


> 


0.000044 


101 


0.000275 


> 


0.000064 


1009 


0.002636 


> 


0.000156 


10007 


0.088889 


> 


0.002733 


20011 


0.183722 


< 


0.337514 


30011 


0.294120 


< 


0.816416 


100003 


1.011050 


< 


10.477100 


200003 


2.117640 


< 


49.372000 


300007 


3.574630 


< 


121.383000 


1000003 


12.647500 


< 


1373.750000 



Table 1: Time each of two programs takes to compute whether p is a Lerch prime. 



2.4 Generalizations 

Euler and Gauss extended Fermat's little theorem and Wilson's theorem, respec- 
tively, to congruences with a composite modulus « — see [ 16, Theorems 7 1 and 1 29] . 
The corresponding generalizations of Fermat and Wilson quotients and Wilson 
primes are called Eider quotients q n (a), generalized Wilson quotients w„, and Wil- 
son numbers n \ w„ (see [28, sequences A157249 and A157250]). (The w„ are not 
called "Gauss quotients;" that term appears in the theory of hypergeometric func- 
tions.) In 1998 Agoh, Dilcher, and Skula [1, Proposition 2.1] (see also Dobson [8] 
and Cosgrave and Dilcher [3]) extended Lerch's formula to a congruence between 
the q„(a) and w„. 

Armed with these facts, one can define generalized Lerch quotients l n and Lerch 
numbers n \ £,,. But that's another story for another time. 



2.5 Open Problems 

To conclude this section, we pose some open problems. 

1. Is £5 = 13 the only prime Lerch quotient? 

2. Is there a fifth Lerch prime? Are there infinitely many? 

Of the 78498 primes p < 10 6 , only four are Lerch primes. Thus the answer to the 
next question is clearly yes; the only thing lacking is a proof! 
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3. Do infinitely many non-Lerch primes exist? 

As the known Lerch primes 3, 103,839,2237 are distinct from the known Wilson 
primes 5, 13,563, we may ask: 

4. Is it possible for a number to be a Lerch prime and a Wilson prime simultane- 
ously? 

Denoting the nth prime by p n , the known Wilson primes are P3,P6,Pi03- The 
primes among the indices 3,6, 103, namely, 3 and 103, are Lerch primes. This leads 
to the question: 

5. If p„ is a Wilson prime and n is prime, must n be a Lerch prime? 

The answer to the converse question — if n is a Lerch prime, must p n be a Wilson 
prime? — is no: p^ and P2231 lie strictly between 563 and 5 x 10 s , where according 
to [4] there are no Wilson primes. 

In connection with Problem 5, compare Davis's "Are there coincidences in math- 
ematics?" [5] and Guy's "The strong law of small numbers" [14]. 



3 Fermat- Wilson quotients and the WW primes 2, 3, 14771 

Suppose that a prime p is not a Wilson prime, so that p does not divide its Wilson 
quotient w p . Then in the Fermat quotient q p (a) of p base a, we may take a = w p . 

Definition 3. If p is a non- Wilson prime, then the Fermat-Wilson quotient of p is 
the integer 



For short we write 



w' - 1 



The first five non- Wilson primes are 2,3,7, 1 1, 17. Since W2 = W3 = 1, wj = 103, 
and w\\ = 329891, the first four Fermat-Wilson quotients are g2 = g3 = 0, 

103 6 - 1 

gl = = 170578899504, 

and 

329891 10 -1 



'11 = 



11 

= 1387752405580695978098914368989316131852701063520729400 



[28, Sequence A197633]. The fifth one, gn, is a 193-digit number. 



Lerch Quotients and Primes, Fermat-Wilson quotients, and the WW Primes 2, 3, 14771 9 

3.1 The GCD of all Fermat-Wilson quotients 

We saw that at least one Lerch quotient and seven Wilson quotients are prime num- 
bers. What about Fermat-Wilson quotients? 

Theorem 2. The greatest common divisor of all Fermat-Wilson quotients is 24. In 
particular, q p (w p ) is never prime. 

Proof. The prime factorizations of q p (w p ) = g p for p = 7 and 1 1 are 
g 7 = 2 4 -3 2 -13-17-19-79-3571 

and 

gu = 2 3 -3-5 2 -37-61 -71 -271 -743 ■ 2999 ■ 89671 -44876831 
• 743417279981-7989680529881. 

Since g2 = gi = 0, we thus have 

GCD(g 2 ,g3 1 g7 1 gii) = 2 3 -3 =24. 

To complete the proof, we show that 24 divides g p whenever p > 3. Since 

pw p = (p-l)\ + l, 

it is clear that if p > 5, then pw p , and hence w p , is not divisible by 2 or 3. As even 
powers of such numbers are = 1 (mod 8) and = 1 (mod 3), and so = 1 (mod 24), 
it follows that pg p (= Wp 1 — 1), and hence g p , is divisible by 24. □ 



3.2 Wieferich primes base a 

Given an integer a, a prime p is called a Wieferich prime base a if the supercongru- 
ence 

a p - 1 = 1 (mod p 2 ) (12) 
holds. For instance, 1 1 is a Wieferich prime base 3, because 

3 10 - 1 =59048= ll 2 -488. 

Paraphrasing Ribenboim [25, p. 264], it should be noted that, contrary to the 
congruence a p ~ l = 1 (mod p) which is satisfied by every prime p not dividing a, 
the Wieferich supercongruence (12) is very rarely satisfied. 

When it is, p cannot divide a, and so the Fermat quotient q p (a) is an integer. In 
fact, (1) shows that a prime p is a Wieferich prime base a if and only if p does not 
divide a but does divide q p {a). 
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In 1909, while still a graduate student at the University of Miinster in Germany, 
Wieferich created a sensation with a result related to Fermat's Last Theorem: If 
X p + y p = z p , where p is an odd prime not dividing any of the integers x,y, or z, 
then p is a Wieferich prime base 2. One year later, Mirimanoff proved that p is also 
a Wieferich prime base 3. (See [6, pp. 1 10-111], [26, Chapter 8], and [30, p. 163].) 

The only known Wieferich primes base 2 (also simply called Wieferich primes) 
are 1093 and 3511, discovered by Meissner in 1913 and Beeger in 1922, respec- 
tively. In 2011 Dorais and Klyve [9] computed that there are no others up to 
6.7 x 10 15 . It is unknown whether infinitely many exist. (Neither is it known whether 
there are infinitely many non-Wieferich primes base 2. However, Silverman has 
proved it assuming the afoc-conjecture — see his pleasantly-written paper [27] .) Like- 
wise, only two Wieferich primes base 3 (also known as Mirimanoff primes) have 
been found, namely, 11 and 1006003. The second one was uncovered by Kloss 
in 1965. An unanswered question is whether it is possible for a number to be a 
Wieferich prime base 2 and base 3 simultaneously. (See [15, section A3] and [25, 
pp. 263-276,333-334].) 

For tables of all Wieferich primes p base a with 2 < p < 2 32 and 2 < a < 99, see 
Montgomery [22]. 



3.3 The Wieferich-non- Wilson primes 2, 3, 14771 

Let us consider Wieferich primes p base a where a is the Wilson quotient of p. 

Definition 4. Let p be a non- Wilson prime, so that its Fermat- Wilson quotient 
q p {w p ) is an integer. If p divides q p (w p ) — equivalently, if the supercongruence 

w p - x = 1 (mod p 2 ) (13) 

holds — then p is a Wieferich prime base w p , by definition (12). In that case, we 
call p a Wiefe rich-non-Wilson prune, or WW prime for short. 

For the non- Wilson primes p = 2,3,7,11, 17,19,23,29,31,37,41,43,47,53,59, 
61,67,71,73,79,83,..., the Fermat- Wilson quotients q P (w p ) =gp are congruent 
modulo p to 

g p = 0,0,6,7, 9, 7, 1,6,18, 17, 30,11, 25, 30,24, 46,64,16,18, 4,29,... (mod p) 

[28, Sequence A197634]. In particular, 2 and 3 are WW primes. But they are triv- 
ially so, because g2 and g^ are equal to zero. 

Is there a "non-trivial" WW prime? Perhaps surprisingly, the answer is yes but 
the smallest one is 14771. In the next subsection, we give some details on using a 
computer to show that 14771 is a WW prime. It is "non-trivial" because g 14771 ^ 0. 
In fact, taking logarithms, one finds that 
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11 



g 14771 — 



( 



14770! + ! \ 14770 
14771 J 



1 



> 10 



8x10 



14771 



so that the number g 14771 has more than 800 million decimal digits. 

3.4 Computer search 

To search for WW primes, one can use a computer to calculate whether or not a 
given prime p satisfies condition (13). Explicitly, if the number 



is equal to 1, then p is a WW prime. 

Mathematical function Mod[a,m] can compute (14) when p is small. But if p 
is large, an "Overflow" message results. However, it is easy to see that in (14) one 
may replace (p— 1)! with (p— 1)! mod /? 3 , a much smaller number. 

For example, it takes just a few seconds for a program using the code 

If [PowerModl (Mod[ (p-1) ! , p"3] + l)/p, p-1, p"2] == 1, Print [p] ] 

to test the first 2000 primes and print the WW primes 2,3, 14771 (see [28, Sequence 
A197635]). 

Michael Mossinghoff, employing the GMP library [12], has computed that there 
are no other WW primes up to 10 . 



3.5 More open problems 

We conclude with three more open problems. 

6. Can one prove that 14771 is a WW prime (i.e., that 14771 divides g 14771) without 
using a computer? 

Such a proof would be analogous to those given by Landau and Beeger that 1093 
and 3511, respectively, are Wieferich primes base 2 . (See Theorem 9 1 and the notes 
on Chapter VI in [16], and "History and search status" in [31].) However, proofs 
for Wieferich primes are comparatively easy, because (high) powers are easy to 
calculate in modular arithmetic, whereas factorials are unlikely to be calculable in 
logarithmic time. 

7. Is there a fourth WW prime? Are there infinitely many? 

Comments similar to those preceding Problem 3 also apply to the next question. 




(14) 



8. Do infinitely many non-WW primes exist? 
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Is it possible to solve Problem 3 or Problem 8 assuming the afrc-conjecture? (See 
the remark in Section 3.2 about Silverman's proof.) 
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